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SYMBOLS 


a,b,c components of in the x, y, and z-directions 


[B] , [Ba ] , [Bg] transformation matrices describing the bending at P for the bend- 
twist sequence with no pretwist, for the twist-bend sequence with no 
pretwist, and for the pretwist-bend-twist sequence, respectively 


n R 

n t~ n n” n b 

n x‘ n y’ n z 

^B^yB^zB 


axis for the bending rotation 

tangential, normal, and binormal unit vectors to a space curve 
unit vectors in the x, y, and z-directions 

intermediate triad of unit vectors corresponding to n x -ny-n z after 
the pretwist and bend rotations of the pretwist-bend-twist sequence 


^(f) ^y^ 


intermediate triad of unit vectors corresponding to n x -ny-n z after 
the twist rotation in the twist-bend sequence 


n x3~ n y 3 n z3 

n £“ n rT n C 

P 


triad of mutually perpendicular unit vectors corresponding to 
n x -ny-n z after deformation 

triad of unit vectors corresponding to n x -ny-n z after pretwist 
arbitrary point on the elastic axis before deformation 


P * 


point P after deformation 


R 


vector position of P 


R* 


vector position of P* 


AR 


vector difference between R and R* 


[T] 

W] ,Wq} ,Wa) 


transformation matrix between n c -n -n, and n_, -n. -n_, 

c, T| s X3 y 3 ^ 3 

transformation matrix between n x -ny-n z and n x3 -ny3-n 23 for no pre- 
twist, initial pretwist, and combined twist cases 


[^ B ] » 


transformation matrix between n -n -n z andn x3 -n -n 23 determined 
with the two-term sequence for t£e lero and nonzero-pretwist cases 


u , v , w displacements at P in the x, y, and z-directions 

x-y-z axes characterizing blade geometry before pretwist and deformation 

where x is directed along the elastic axis 


x ,-y,- 2 , 


curvilinear coordinates of deformed blade where x 3 is directed along 
the elastic axis 


a 


sum of 6 and <f> a 



3*3Q»3 a flap rotation about -y or updated -y-axis for the no pretwist, for 

the initial pretwist, and for the combined twist cases, respectively 

e some small nondimensional value corresponding in magnitude to the slope 

of the deformed blade 


S-n-C 


curvilinear coordinates of a pretwisted blade with £ directed along 
the elastic axis 


C.Ce.Ca 


e 

K 

T 

$ 


*e 

5 


lag rotation about the z-axis or updated z-axis for the no pretwist, 
for the initial pretwist, and for the combined twist cases, 
respectively 

pretwist angle 

bending curvature of a space curve 
torsional curvature for a space curve 
bending rotation angle 

deformation twist variable defined in terms of Euler-like sequence for 
the no pretwist, for the initial pretwist, and for the combined twist 
cases, respectively 

elastic twist 

deformation twist variable defined in terms of the two-term sequence 
for the zero- and nonzero-pretwist cases 


co x3 ,a)y 3 ,o) 23 components of u> in the x 3 , y 3 , and z 3 -directions 
0 ) curvature vector of the blade 


a 


angular velocity of rotor 


a 


curvature vector of space curve 


Subscripts : 

1,2, 3, 4, 5, 6 refer to the flap-lag-twist (1), flap-twist-lag (2), twist-flap-lag (3), 
lag-flap-twist (4), lag-twist-flap (5), and twist-lag-flap sequences (6) 


iv 



SUMMARY 


Adequate modeling of elastic helicopter rotor blades experiencing moderately 
large deformations requires a nonlinear analysis. This analysis must be based on an 
appropriate description of the blade’s deformation geometry, including elastic bending 
and twist. Built-in pretwist angles complicate the deformation process and its defi- 
nition. In this study, relationships between the twist variables associated with 
different rotation sequences are listed, as well as corresponding forms of the trans- 
formation matrix. Included are relationships between the twist variables associated 
with, first, the pretwist applied initially and, second, the pretwist combined with 
the deformation twist. Many of the corresponding forms of the transformation matrix 
for the two cases are also listed. Moreover, it is shown that twist variables con- 
nected with the combined-twist treatment can be related to those in which the pretwist 
is applied initially. A method is outlined for determining these relationships, and 
some results are given. Additionally, a procedure is demonstrated for evaluating the 
transformation matrix that eliminates the Euler-like sequence altogether. The result- 
ing form of the transformation matrix is unaffected by rotation sequence or pretwist 
treatment. 


INTRODUCTION 


To derive the equations of motion for an elastic rotor blade, one usually begins 
by relating the blade’s instantaneous deformed geometry to that of the undeformed 
blade. It is common practice in the helicopter literature to evaluate the transfor- 
mation matrix between the deformed and undeformed states using an , Euler-like sequence 
of three successive rotations (e.g., refs. 1-4). Furthermore, in the helicopter 
literature the built-in twist of the undeformed blade is often included with the 
deformation twist in the rotation sequence (e.g., refs. 1 and 2). For linear mathe- 
matical models (ref. 5), the order of rotation and the handling of the pretwist do 
not affect the final form of the transformation matrix. However, in nonlinear analy- 
ses, the final form of the transformation matrix — and subsequently the derived equa- 
tions of motion — will vary depending on the deformation sequence and pretwist 
treatment. 

Adequate modeling of elastic rotor blades experiencing moderately large deforma- 
tions requires a nonlinear analysis. In basing that analysis on an appropriate 
description of the deformation geometry, the question arises of what form to use for 
the transformation matrix and how it relates to other forms. Since the built-in pre- 
twist angle is physically present before the deformation, applying the pretwist 
before the blade bending and twisting deformations is unquestionably correct. How- 
ever, combining the pretwist with the deformation twist has not been shown to be less 
correct. This combination of twists usually simplifies the derivation and resulting 
expressions — but the validity of the approach should be examined. 

Kaza and Kvaternik (refs. 2, 6, and 7) developed the curvatures, transformation 
matrices, and equations of motion for a rotating elastic blade^ using the flap-lag-twist 


and lag-flap-twist rotation sequences concurrently. The authors assumed that the 
twist variables used in the two sequences were identical and equal to the elastic 
twist. Note that the elastic twist is that part of the deformation twist associated 
with the torsion of the blade. Bending can also contribute to the deformation twist. 
In not distinguishing between the elastic and deformation twists, Kaza and Kvaternik 
concluded that the two different sets of equations were not equivalent. The authors 
suggested interpreting 0 the results as M two different nonlinear approximations of a 
given physical system” (ref. 2, p. 81). 

Rosen and Friedmann (ref. 8, p. 163) concluded that the orientation of the triad 
of unit vectors characterizing the deformed blade can differ by second-order terms 
that depend on the rotation sequence. Therefore, these authors chose one sequence 
and were consistent throughout the derivation (ref. 3, p. 45 or ref. 8, p. 163). 

Petersen (ref. 9) showed that the change in orientation of a rod with zero 
elastic twist could not be properly described by either a flap-lag or a lag-flap 
rotation sequence. Since zero elastic twist does not mean zero deformation twist, 
it is not surprising that the flap-lag and lag-flap sequences are inadequate. The 
resulting orientations differed from each other and from the actual orientation by 
rotations about the deformed elastic axis. The actual orientation could be deter- 
mined by either a flap-lag-twist or a lag-flap-twist sequence in which two different 
twist variables were evaluated with the zero elastic twist condition. Petersen 
extended the results to include elastic and built-in twist. 

Hodges et al. (ref. 10) showed that by regarding the twist variable associated 
with each Euler-like sequence as different, it was possible to find relationships 
between the twist variables for the different sequences. Using these relationships, 
it was then possible with the appropriate change of variable to make any one form of 
the transformation matrix identical to any other. The authors reasoned that regard- 
less of the rotation sequence used in the analysis, an elastic beam with a given set 
of end conditions and loads is a single physical system. It was argued that the 
transformation matrix and curvatures were unique, but that the form varied according 
to the choice of variables — specifically, the choice of the twist variable. Two 
methods for determining the relationships between the different twist variables were 
employed. The principal approach used in reference 10 involved equating expressions 
for an integrated form of the torsional curvature. Another was also outlined 
(ref. 10, pp. 42 and 46) — to equate components of the various forms of the trans- 
formation matrix. Hodges et al. determined a relationship between the twist variables 
associated with the flap-lag-twist and the lag-flap-twist rotation sequences. 

Several authors have made it a point to apply the pretwist before the deformation 
rotations (refs. 3 and 4). However, the arguments presented seem to rest more on 
"realism” than on the inappropriateness of combining the twists (ref. 4, pp. 40 
and 41). In response, Hodges et al . (ref. 10, p. 20) claimed that the pretwist could 
be combined with the deformation twist in the flap-lag-twist and lag-flap-twist, but 
not in the flap-twist-lag, twist-flap-lag, lag-twist-f lap , or twist-lag-flap rotation 
sequences. The authors argued that the pretwist, elastic twist, and the deformation 
twists are all rotations about the same axis for the flap-lag-twist and lag-flap- 
twist sequences. 

The results of reference 10 for the case of zero pretwist are generalized and 
extended in this paper. Relationships between the twist variables associated with 
different rotation sequences are listed, as well as corresponding forms of the trans- 
formation matrix. Included are relationships between the twist variables associated 
with, first, the pretwist applied initially and, second, the pretwist combined with 
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the deformation twist. Many of the corresponding forms of the transformation 
matrix for the two cases are also listed. Moreover, it is shown that twist variables 
connected with the combined-twist treatment can be related to those in which the pre- 
twist is applied initially. A method is outlined for determining these relation- 
ships and some results are given. Additionally, a procedure is demonstrated for 
evaluating the transformation matrix that eliminates the Euler-like sequence alto- 
gether. The resulting form of the transformation matrix is unaffected by rotation 
sequence or pretwist treatment. 

The emphasis in this paper is on generalizing some of the arguments presented in 
reference 10, validating the approach of lumping the pretwist with the deformation 
twist, offering an alternative approach for evaluating the transformation matrix, and 
providing the means to compare sets of equations that have been derived using differ- 
ent rotation sequences and different pretwist treatments. 


DEFORMATION GEOMETRY 


Assume that the elastic axis of the undeformed blade lies along the x-axis. 
Now, consider the blade frozen in some deformed state. At that instant, the elastic 
axis lies along a unique curve in three-dimensional space (fig. 1).. Before deforma- 
tion, the position of a point, P, on the elastic axis is 


R = xn x 


After deformation, the position of P, now P*, is 

R* = R + AR 

= (x + u)n x + viiy + wn z 

At any given instant, the displacement variables, u-v-w, are functions of position 
along the undeformed elastic axis, x. As a result, 

dR* = dx [ (1 4- u’)n x + v'n y + w'n z 

where ( ) 1 indicates differentiation with respect to x. 

The curvilinear coordinate, x 3 , is measured along the deformed elastic axis; 
is the unit vector tangential to the deformed elastic axis. Therefore, 

dR* = dx n 

3 x 3 

The magnitude of dR* is equal to dx 3 and also to 

II dR* II = (dR* • dR*) l/2 

= dx[(l + u') 2 + (v') 2 + (w') 2 ] 1 / 2 

Equating expressions for II dR* li yields a relationship between dx and dx 3 : 
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dx 


= [d + u ') 2 + (v ') 2 + (w’) 2 r l/2 

3 


Then, n X3 can be expressed as 

"xa = ^7 [(1 + u ’ )n x + v,n y + w ' n zl 


Let 


n x3 = an x + bn y + cn z 

where a, b, and c represent the previously derived components of n^. The relation 

a 2 + b 2 + c 2 = 1 

will be used repeatedly throughout the paper. 

In helicopter-blade aeroelastic analyses, it is usually assumed that although 
deflections v and w can be moderately large, slopes v 1 and w’ are small — of the 
order of some e in size. In general, since rotor blades are long and slender, it 
is assumed that u f is smaller — of the order of e 2 . Using these guidelines and 
retaining terms through 0(e 2 ), 

a s 1 - j (v’) 2 - 

b s v’ 

C 2 W 1 

Thus, a is 0(1), and b and c are 0(e). To 


j (w’) 2 


0(e 2 ), a can also be expressed as 


TRANSFORMATION MATRIX 


In general, fibers at P along the x-, y-, and z-directions do not remain 
mutually perpendicular, but distort in proportion to the shear strains (ref. 11, 
p. 23). Shear deformations are neglected in helicopter-blade analyses, however, and 
the local n x -Uy-n z triad is assumed to rotate into another triad of mutually per- 
pendicular unit vectors, n x3 ™ n y 3 “n Z3 . Note that if shear effects are included, the 
change in orientation at P can be broken into a rigid body part plus a distortion 
part due to the shear effects. The rigid-body rotation is described in terms of a 
transformation matrix, , where 
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Three of the nine components of \<y\ are known; the other six must be determined. 

The variables that describe the displacement at P specify the direction of 
n x3 and the plane in which n 3 and n Z3 fall. If the first row of [&] is as shown, 
the final orientation described by any orthonormal [ff’] will differ at most by a rota- 
tion about n X3 . Therefore, whatever the scheme used to evaluate the remaining com- 
ponents of [ST] , if orthonormality is maintained, \ST\ will describe the desired 
orientation to within some twist angle. 

For the analytical model to include the effects of elastic torsional deformation, 
an additional degree of freedom is required. Consider incorporating a rotational 
degree of freedom about the elastic axis — one that is defined in the process of 
evaluating the transformation matrix. This deformation-twist variable includes the 
elastic twist along with some part of the bending twist. The part of the bending 
twist that is included depends on the sequence used to evaluate the transformation 
matrix. By including a deformation-twist variable, the analytical model is given 
sufficient flexibility to properly orient n 3 and n Z3 . Twists like those developed 
in reference 9 to correct the lag-flap and flap-lag results for a rod with zero elas- 
tic twist are simply lumped into the deformation-twist variables. For the six Euler- 
like sequences examined here and in reference 10, along with the alternative approach 
developed here, the different twist variables and associated forms of the transfor- 
mation matrix should offer equivalent descriptions of the deformed blade. In this 
context, the reasoning of Hodges et al. (ref. 10) becomes appropriate. 

Whether or not the blade has built-in twist, the orientation of n x3 -ny 3 -n Z3 
is still described in terms of i^-ny-iiz by an orthonormal matrix with a, b, and c 
in the top row. Although this transformation matrix includes a pretwist contribution, 
the conclusions of the last two paragraphs remain valid. 


EVALUATING TRANSFORMATION MATRICES AND RELATING TWIST VARIABLES 

No Pretwist 

The lag-flap-twist sequence corresponds to the series of rotations (fig. 2) : 

1. about the z-axis 

2. about the updated y-axis 

3. about the twice updated x-axis 

The transformation matrix associated with the lag-flap-twist sequence is 
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w = 


cos 3 ^ cos £ ^ cos 6 4 sin? 4 sing^ 

-cosd), sin?, - sind>, sing, cos?, cos<j>, cos?, - sin<i> sing, sin? sinf cosg 

T 4 4 4 4 4 4 4 4 4 4 4 4 

sin<}> 4 sin^ 4 - cos(ji 4 sing 4 cos£ 4 -sin<t> 4 cos£ 4 - costj) 4 sin$ 4 sin£ 4 cos(j) 4 cosg 4 

f 

By equating the components of the top row to a, b, and c, respectively, expressions 
for the sines and cosines of and g 4 can be determined. These can then be sub- 
stituted back into ] . By defining the twist variable in terms of (j> 4 , ] is 

completely determined. An alternative approach for defining the twist variable is 
discussed later in this report. 


Based on the rotation scheme discussed above, there are six possible sequences: 

1. Flap-lag-twist 

2. Flap-twist-lag (P 2 " ( I> 2 “C 2 ) 

3. Twist-flap-lag (<f> 3 -p 3 -£ 3 ) 

4. Lag-flap-twist (C 4 “S 4 -<}> 4 ) 

5. Lag-twist-flap (£ 5 “<f> 5 “B 5 ) 

6. Twist-lag-flap (<f> 6 ~£ 6 -3 6 ) 

In terms of a, b, c, and the appropriate twist variable, the final form of [ST] asso- 
ciated with each rotation sequence is listed in table 1. 


Rotor blades are relatively stiff in torsion, and twist deformations are usually 
assumed to be small — of 0(e 2 ) in size. Retaining all terms through 0(e 2 ), there 


are two distinct 

forms for . 

The flap-lag 

form 

can be 

expressed 

as 


"l - 

b 2 /2 

- c 2 /2 

b 


c 

- 



[^.1 = 


-b - 

1 

- b 2 /2 - <J>?/2 


4>i - be 

» 

i = 1,2,3 


- 

-c + 

b<j>i 


1 - 

c 2 /2 - 

^±/ 2 J 



and the lag-flap 

form 

as 








1 - 

b 2 /2 

- c 2 /2 

b 


c 

- 





-b - 

c *3 1 

- b 2 /2 - <J)j/2 


*3 


9 

j = 4,5,6 



-c 4* 

b *j 

_ +j ‘ bc 

1 - 

c 2 /2 - 

*]n_ 




For the linear case, in which terms smaller than 0(e) are neglected, all six trans- 
formations reduce to the single form. 
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1 b c 

-b 1 4> k 

L-c -<t> k 1 , 


k = 1 , 2 , 3 , 4 , 5 , 6 


Relationships between the six twist variables introduced above are listed in 
table 2 . Included are relationships of each of the other five twist variables to 
4 > 1 . Manipulating these expressions, one can also find relationships between the five 
The relationships of the lag-flap-twist variables, <|> 5 and <j> 6 , to <j> 4 are listed. 
Several others are added for use later. These expressions were derived by equating 
components of the different forms of W\ . It is easily verified that by substituting 
the expression for one twist variable in terms of another, identical transformation 
matrices are produced. 


Initial Twist 

The pretwist, 0 , is built-in about the blade f s elastic axis before deformation. 
This pretwist can be a function of x but not of time. Figure 3 illustrates the 
undeformed geometry of the blade. The tangent to the elastic axis after deformation 
can still be described as 


"xa = an x + bn y + cn z 


Relating the deformed geometry to that of the undeformed blade 


m 


X3 




n. 


y3 

K. n zJ 


= [T]. 


rn f| 

"n > = 




■A 

L 


L L 3 



where 


A = a 

B = b cos 0 + c sin 0 
C = -b sin 0 + c cos 0 


The lag-flap-twist (£ q 4 — ^ ^ g 4 ) sequence is illustrated in figure 4 . To evalu- 
ate [T 4 ], substitute A for a, B for b, C for c and 4>0 4 for <£ 4 in [^ ] . 

In the same manner, the forms of [T] associated with the flap-lag-twist (6 Ql -c ei -d> 0 ) 
flap-twist-lag » twist-flap-lag (4103-303-403), lag-twist-flap 1 


The triad, n x3 -n y3 -n 23 . 

can also be 

related 

to 

tly- 

-n z by 


rn x 3 * > 


rn x^ 


"1 

0 

0 “ 

< 

“ya 

>• - 

n y 

> = [T] 

0 

cos 0 

sin 0 


^ n Z 3 J 




0 

-sin 0 

cos 


rn x 

n y 


Several forms of W§] are listed in table 3 . 
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Again, when retaining terms to 0(e 2 ), the six forms of Wq] reduce to two. 
flap-lag form is expressed as 


The 


[^ 0i ] 


r 2 2 

1 ^ — 

2 2 


-b cos 0 - c sin 0 -(<f>g^ + bc)sin 0 


(<t>Qi - bc)cos0 


+ 4 > 0 j[(b sin 0 - c cos 0 ) + ll-b 2 +y(b cos 0 + c sin 0) 2 + l-c 2 +y(b cos 0 + c sin 0 ) 


[i- b2+ i 

% 


[ lc2+ i 

‘] 


b sin 0 - c cos 0 -< ^ 0 i cos 9 


-4>0i sin© 


+ 4 > g i (b cos 0 + c sin 0 ) 


- 1 1 — -^-(b sin 0 -c cos 0) 2 


t 

+Sil . e 

— sm 6 


+ l--^-(b sin 0 - c cos 0 )' 


H 


i = 1,2,3 


and the lag-flap form as 




u 2 2 

b c 


1 2 2 

b 

c 

-b cos 0 - c sin 0 

- 4 > e j sin 0 

(f)0j COS0 

+ (fjQj (b sin 0 - c cos 0 ) 

+ |l-y(b cos 0 +c sin 0) 2 

+ ll-y(b cos 0 +c sin 0) 2 


♦y , 

2 cos 0 

L 

- 2 J jsi„e 

b sin 0 - c cos 0 

+ bc)cos 0 

-( 4 » 0 j - bc)sin 0 


+ <}) 0 j (b cos 0 +c sin 0 ) 




«t>L 


- |l-b 2 +y(b sin 0 - c cos 0) 2 + 


[ 1-c2+ i 


-^-(b sin 0 - c cos 0 ) ^ 


- sin 0 

A, 2 T 

- , 


2 


cos 0 


j = 4,5,6 
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Neglecting terms smaller than 0(e), all six forms of Wq] reduce to 


1 


b 


- 


-b cos 0 - c sin0 cos0 - <{)q^ sin0 

b sin0 - c cos0 -sin0 - cos0 


c 

sin0 + 4 >q^ cos0 ; 
cos0 - <f>ek sin0 


k = 1,2, 3, 4,5, 6 


Note that no assumption has been made about the size of 0. 

Relationships between the six twist variables introduced above can be determined 
by equating components of the different forms of [T] . Equivalently, the expressions 
derived for the case of no-pretwist (table 2) can be extended to this case by substi- 
tuting A for a, B for b, C for c, and <f>Q for <J>. Some revised expressions 
are listed in table 4. The relationships now include a 0-dependence. 


Combined Twists 

Let the transformation matrix, [F a ] » describe the difference in orientation 
between n x " n y“ n z ancl n x 3 ~ n y 3 “ n Z 3 * As before, the components of the top row are 
equated to a, b, and c, respectively. The pretwist is combined with the deformation 
twist in the rotation sequence used to evaluate the rest of [^J. Six rotation 
sequences are considered: 


1. 

Flap-lag-twist 


2. 

Flap-twist-lag 

(8 a2 ~ a 2 - ^a2^ 

3. 

Twist-flap-lag 

( a 3 - ^a3 - ^a3^ 

4. 

Lag-f lap-twist 

(4 a i* - ^a4 -a 4^ 

5. 

Lag-twist-flap 

^a5 -a 5 - ®as^ 

6. 

Twist -^lag-flap 

( a 6~4a6 - ^a6 ) 


where 

a i = 0 + ^ai ; 1 = 1» 2 ,3,4,5,6 


A different form of [^] is associated with each rotation sequence. The results for 
[STj] given in table 1 can be used for [^ a ±] by substituting a-^ for (j)^. 


In table 5, expressions for are listed in which all terms greater than 

0(e 2 ) are neglected. Neglecting terms smaller than 0(e) reduces the six forms of 
[3T a ] to 





1 


b 



c 



W«i1 - 

-b 

COS0 

- c sin0 

cos 0 

~ *ai 

sin0 

sin0 

+ <t> . 
v ai 

cos0 

> 


_ b 

sin0 

- C COS0 

-sin0 

" *oi 

COS0 

COS0 

^ai 

sin0 



i = 1,2, 3, 4, 5, 6 
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Note that the linear expression above is identical to that derived for the case of 
initially applied pretwist. 


By equating components of the different [<^J, relationships can be determined 
between the different a-variables. Since a is not necessarily of 0(e), the 
results to 0(e 2 ) that are listed in table 2 cannot be extended to the present case. 
However, the other results in table 2 can be extended by replacing <J>^ with a-^. 
Some modified expressions are listed in table 6. As before, the different forms of 
[S a ] can be reduced to one by the appropriate change of the twist variable, a. 


The trigonometric identity, 

tan(<J> a + 6) 


tan<j> a + tan0 
1 - tan<j> a tan0 


can be used to relate the <f> a -variables directly. Some expressions are listed in 
table 7. Although the relationships are largely 0-dependent, the relationship 
between 4> al and <j> a4 is identical to the zero pretwist case. Recall that. <J> ai is 
associated with the flap-lag-twist sequence and <J> ^ with the lag-flap-twist 
sequence. This result is a direct consequence of the matrix identity 


*1 

0 

0 


“i 

0 

o - 


“i 

0 

0 

0 

cosa 

sina 

= 

0 

COS0 

sind 


0 

cos<f> 

sin<£ 

.0 

-sina 

cosa_ 


_0 

-sin0 

cosB 


0 

-sincj) 

COS(J)_ 


Both [< 5^1 and [S' 4 ] can be expressed as products of the 0-matrix and matrices inde- 
pendent of 0. Therefore, in equating [S~ al ] and the expressions are indepen- 

dent of 0. 


RELATING TWIST VARIABLES ASSOCIATED WITH THE INITIAL PRETWIST 
AND COMBINED PRETWIST TREATMENTS 


Relationships between twist variables have been determined for the case in which 
the pretwist is applied initially and for the case in which the pretwist is combined 
with the deformation twist. It is possible, by equating components of the appropriate 
forms of the transformation matrix, to relate twist variables associated with one pre- 
twist treatment to those associated with the other. Equivalently, results derived in 
previous sections can be used to determine the desired relationships. 


For the twist-flap-lag and twist-lag-flap rotation sequences, lumping the pre- 
twist with twist in the sequence does not alter the final form of the transformation 
matrix; that is, 


[S’ A and [S' 1 are identical in form to [SI ] and 1, respec- 

r a3 J t a6 0 3 J A c J * r 

t ively . Therefore , 


06 - 


^a3 ^03 

^a6 “ ^06 
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Consider one of the six forms of [^ a ] — say The components of [<^ ai ] are 

functions of a, b, c, and c^. This form of the transformation matrix can be made 
identical to [^, 3 ] by substituting for a 3 in terms of a 3 where (from table 6) 

a tana 3 

tana. = ; 

1 (1 - b 2 ) - be tana 3 

By definition, 

°3 = 6 + *03 

and, since [^* 3 ] is identical to [^ 03 ], 

^063 ~ ^03 

Substituting for a 3 and using trigonometric identities, [<^t 3 ] can be expressed in 
terms of the sines and cosines of <j> a3 and, therefore, of <f>e 3 - Then, using another 
change of variable based on the relationship (from table 4 ), 

[1 - (b cos0 + c sin0) 2 ] tan<)>0 1 
tan<f> 0 3 a _j_ s i n 0 + c C os0)(b cos 0 + c sin0)tan0 1 


what was originally [& al ] can be made identical to [^ x ] * 

Figure 5 illustrates the process with a flow chart. Although the procedure out- 
lined above is indirect, it shows that, in general, any [«^J can be made identical to 
any other ] by a series of previously derived changes of twist variable. The 
sequence of variable changes has been condensed for the two examples given in table 8. 


ALTERNATIVE METHOD FOR EVALUATING THE TRANSFORMATION MATRIX 


In general, a sequence of three Euler-like rotations is used to evaluate six 
unknown components of the transformation matrix. Suppose the rotation sequence was 
defined in terms of two sequential rotations instead of three. The overall rotation 
could be broken into a bending part and a twisting part. As before, 

°X3 = aQ x + bn y + cn z 

where a, b, and c are functions of the displacement at P and a 2 + b 2 + c 2 = 1 . 
To characterize the bending rotation, an axis of rotation and an angle of rotation 
about this axis are required. The axis of rotation, n^, is to be perpendicular to 
both n x and n x3 (fig. 6). To determine n^, then. 
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n x x n X3 


nR l|n x * n x3 


. / 2 ”y 


+ 


<b 2 + c 2 ) 1 ' 2 " z 


(b 2 + c 2 ) 1 

+ ^y n y + ^z n z 

The sine and cosine of the angle of rotation, $, are 

cos = n x • n x3 = a 

sin$ = lln x x n x3 ll = (b 2 4- c 2 ) 1 / 2 

Knowing the axis of rotation, the associated transformation matrix can be evaluated 
with the matrix equation (ref. 12, p. 35) 



_ i 

0 

o' 


"4 

Rjjlly 

*x R z 


’ 0 

R z 

-v 

[B] = cos4> 

0 

1 

0 

4 (1 - cos$) 



^ R Z 

4 sin4 

- R z 

0 

Rx 


.0 

0 

1. 



R z*y 

R l . 


-*y 

-*x 

0 _ 


-b a 4* c 2 /I 4 a -be /I 4 a 

-c -bc/1 4 a a 4 b 2 /l 4 

Consider the bend-twist sequence illustrated in figure 7(a) for the case of no 
pretwist. The transformation matrix between n x 3 ~ n y 3 _n z 3 anc * 


t^ B ] = 


10 O' 

0 cos^-j- sin^-p 
LO -sin^ coscf>,jJ 


[B] 


-b cos(j> T - c sin<j> T ( a + 'T+~a) COS ^T ~ TTa" s ^ n( ^T ( a l^Ta’) sin ^T _ l^a cos % 


j^bsin(f) T -c cos(j> T -|a4 - Y ^ Ta COS( ^T 4 cos^j 4 sintft^j 

Now consider the twist-bend sequence shown in figure 7(b). For this case, the bending 
transformation matrix represents a rotation between n X( j ) -ny ( j ) -n Z( j ) and n X3 -ny 3 -n Z3 . The 
rotation axis can be described in terms of components in the n x<j>“ n y<j)~ n Z(}> directions, 
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[• 


b sinijirj. - c cosij)^ b cos<f>-ji 


(b' 


+ c 2 ) 172 " 


(b 2 + 



The expression for $ is unchanged and the associated transformation matrix becomes 


(BJ = 


-b cosc^-c sin<() T 
b sin^rj, - c cos^ 


a + - 


b cos4>/j.+ c sin<J>rj, 

(b sin<|>j. -c cos^) 2 


-b sin4>£+c cos<}>^ 

(b cos(J)rp+c sin4>.p) (b sin$<p-c cos^^) 


1 + a 

(b cos<j>,j+ c sin<{)^)(b sin^^-c cos^) 
1+ a 


1 + a 


a + - 


(b cos^+c sin^)^ 
1 + a 


Then, the transformation matrix between i^-ny-^ and n x3 -ny 3 ~n Z3 becomes 

noon 


[& B ] - tv 


0 cosfyj sincf)^ 
LO -sincf)r j, cos^J 
a 


c 



/ , c 2 > 

\ r bC 


*b cos(J>r^-c sin<f>j 

( a + 1 + a> 

)cost > T - 1 + a 

S 


/ c 2 ' 

y . , be 


b sin<f)^-c cos(f>rj. 

_ ( a+ l + a> 

|sin<() T 1+a 

COS^fji 


( a+ irir) sln +T-iTi 
( a + ITl) COB +T + iTI 


cos4 

slnfyrp 


It is apparent that there is only one <J>j and associated form for [^3] f° r either 
the bend-twist or the twist-bend sequence. 


Suppose the blade is pretwisted by an amount 0 . For the pretwist-bend-twist 
sequence illustrated in figure 8, the bending transformation matrix represents a 
rotation between n^-n^-n^ and n^-ny^-n^. The rotation axis is 


nr 


-[■ 


b sin0 - 


(b 


2 + c 2 ) 


c cos0 
1J2 


b cos0 


(b 2 + c 2 ) 1 


+ c sinO 

2 U/2 



The expression for <£> is again unaffected and [ Bq ] is the same as [B^] with 0 sub- 
stituted for Then the transformation matrix between n x -ny-n z and n X3 -ny 3 -n Z3 

is 
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1 


0 


0 


0 


0 


[^■gg] = 0 cos(j> T g siting [Bg] 0 cose sin0 


0 -SincflrTnQ COS<j>n 


= l-b cosa-c sina 


0 -sin0 cose 
b 


L , 

c 2 ,> 

i be 

t 1 + 

1 + a> 

icosu— " 

' 1 + a 

f 

.3 _i_ 

c 2 > 

be 

a t 
< 

1 + aJ 

! ® ina - 
' 1 + a 


( a + ITl) sin “-T 


L , / , c 2 \ . be (, o \ , be . | 

I b sina-c cosa -ia + -=— ; — sma-r- — cosa (a + y-; — )cosa + -r— — sma 
L \ 1 + a/ 1 + a \ 1+a/ 1 + a J 

where a = 6 + <J>j 0 . For the pretwist-twist-bend sequence, the transformation matrix 
is identical to that above. From the form of [&■%§], it is apparent that whether the 
transformation is determined by applying 6 first or combined with the twist vari- 
able, the resulting expressions are the same. To 0(e 2 ) in the sense already 
established, 

F i-T-T- b = 1 


t> sin0 - c cos0 


( b 2 fre\ / c 2 *tg\ 

-b cos6 - c sin0 II — ^ ^~lcos6 II — ^ — IsinP 

[^g 0 ] = + 4>x 0 (b sinS “ c cos6) - + -^sinG + ^ TQ - ^=-)cos0 

/ b 2 '^Te\ / c 2 ^T6\ 

d sin 0 - c cos 0 -II — 2 2 ~ j sinQ — 2 2 ~/ COS 

-f (b cos0 + c sin0) - 4- -y^cos0 - (<j> T g - -rpjsin0 

The matrix [<^q] can, with the appropriate change of variable, be made identical 
to any of the previous forms of [<^g] or [^ a ] • The relationship between <|> al and (f>pg 
and between <}> a4 and <£pg are listed here: 

(a + a 2 4- c 2 )tan<(>pg + be 

tan<(> — 

-be tan<})pg + (a 4 a + c ) 

(a + a 2 4- b 2 )tanc()jig - be 

tancf)-.. — — r o o 

be tan<j)pg 4- (a + a 4b) 

To 0(e 2 ), these can be expressed as 


>ai 8 *T0 + T 


*0 4 B ^T0 2 
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The physical significance of cf^ seems clear. It represents the total defor- 
mation twist — the bending, as well as the elastic contribution. Whether this twist 
variable has any real advantages over previously discussed twist variables remains to 
be demonstrated in further studies. 


CURVATURE-DEFINED TWIST VARIABLES 


Another approach for defining the twist variable is through the curvature. As 
before, the transformation matrix is derived in terms of a, b, c, and some 
deformat ion- twist variable. An expression for the torsional curvature is determined 
from the transformation matrix and equated to an elastic-twist rate. Solving for the 
deformation twist in terms of the elastic twist, the expression can be substituted 
back into the transformation matrix. 


To develop the expressions for the curvature, recall the matrix relation 



and note that 



Using a kinetic-energy analogy (ref. 13, pp. 381-385), the curvature vector is defined 
as 


such that 


^ ^xa^s + w y 3 n y 3 + aj z3 n Z3 
dll X3 


dx 3 

= 0) 

x n X3 


W y3 n Z3 

dn y3 

dx 3 

= (a) 

x V 

= -“zs 1 ^ 

+ W X3 n Z3 

dn__ 

Z 3 

dx 3 

0) 

X n z 3 

t>) y3 n yL3 



The components of w can then be solved for 
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dn 


y3 


X 3 


n 


dx 3 “Z3 


-dn . 


X 3 


y 3 


n. 


dx 3 . « 


dn 


X 3 


Z 3 


dx 3 y 3 


and in terms of components in the n x -ny-n z directions, one obtains 

dy 

W X3 = dT z + m 2 m 3 + n 2 n 3 } 


V = (a ' Z 3 + b ’ m 3 + C ’ n 3 ) 


j = 4— (a'£ + b'm + c'n ) 

Z3 dx * 2 2 2 y 


These expressions can also be derived directly from the Serret-Frenet formulas 
for a space curve (ref. 14, pp. 501-503). The Serret-Frenet formulas refer to a 
three-dimensional space curve and the associated tangential (n t ) , normal (i^), and 
binomial (n^) vectors. The derivatives with respect to distance along the space 
curve, s, are 


dn^ 

ds 


= Q x n i 


where ft is the curvature vector. 


i = t ,n ,b 


ft = xn t + 

The unit vector, n X3 corresponds to n t , and ny 3 and n Z3 are in the same plane as 

n 1 n and n^. It is not difficult to relate the two triads in terms of a rotation about 

n t and subsequently derive the expressions for m. Despite statements to the con- 
trary (ref. 7, p. 91), these relations do not rely on the coincidence of n y3 and n Z3 

with n n and — only on the coincidence of with n t . 

For each twist variable and transformation matrix already discussed, there is an 
expression for w. As in the case of the transformation matrix, one form of ca can 
be made like another by appropriately changing the twist variable. The torsional 
curvature, ^ x3 , can be expressed as a function of any one of the previously defined 
twist variables and its derivative, as well as a, b, c, 0 , and their derivatives. 
Ideally, it would be possible to turn the expressions for w x3 around, that is, to 
solve for each of the twist variables in terms of oj , a, b, c, and 0. Then, if 
the resulting expressions were substituted back into the transformations, the differ- 
ent forms of the transformation matrix should reduce to a single form in terms of 
t^x 3 j a, b, c, and 0 . 
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For the case in which the blade has no built-in pretwist, let an elastic-twist 
variable, <j> e , be defined as 


<♦«> ""*! 

In table 9, several relationships are listed between <f> e and the twist variables 
defined through different rotation sequences. Note that in integrating, the deforma- 
tion twist and elastic twist are assumed to be zero at x = 0. For the flap-lag- 
twist and lag-flap-twist sequences, the twist variables, ^ and <j> 4 , respectively, 
can easily be solved for in terms of <J> e and an integral expression of a, b, and c. 
Like 4> x and (J> 4 , the twist variable defined in the previous section, <j)j, can also be 
solved for in terms of <j> e . This is not the case for the other four Euler-like 
sequences. As an example, the relationship between <j > 3 and <j> e is included. 

Finally, when the blade has a built-in pretwist, 0, one way of defining the 
elastic-twist variable, <b e , is 


dfj ( *e + 9 > = w x 3 

In table 10, relationships are listed between <j> e and <Pq i9 <P al > ^ 04 * <f> a 4 > and <j> T Q 
As in the zero-pretwist case, the deformation and elastic twist are assumed to be 
zero at x = 0 . 


COMPARISON OF RESULTS 


Except for reference 10, few results have been published with which the relations 
derived in this paper can be compared. Results listed in reference 10 (pp. 16 and 22) 
are identical to those developed here relating <p al and <|> a4 , <j> al and <p e , and <}> a4 
and <p e . 


CONCLUSIONS 


The deformation geometry of an elastic blade has been investigated. In evaluat- 
ing the transformation matrix between the undeformed-blade triad n x -ny^n z and the 
deformed-blade triad n X3 -ny 3 -n z3 , the direction of n^. and the plane of Uy 3 and 
n Z3 are specified by the deflections of the elastic axis. This is also true when 
there is built-in twist, regardless of whether it is considered to be applied before 
deformation or combined with the deformation twist. Despite different rotation 
sequences and pretwist treatments used to determine the transformation matrix, the 
use of a twist variable defined through the rotation sequence in the mathematical 
model allows each approach to orient n X 3~ n y3~ n z3 identically. 

An alternate approach to using an Euler-like sequence of rotations is discussed 
in which the change in orientation is broken into a bending and a twisting part. The 
results are independent of sequence and pretwist treatment. The twist variable 
defined with this approach represents the total twist deformation. Of the Euler-like 
sequences, the flap-lag-twist and the lag-flap-twist sequences appear to offer the 
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simplest expressions for the transformation matrix and curvature vector. The alter- 
nate bending and twisting approach produces expressions for the transformation matrix 
and curvature vector of about equal simplicity. Whether the alternate approach 
offers any real advantages over the flap-lag-twist and lag-flap-twist sequences 
remains to be demonstrated. 

Relationships between twist variables associated with the six different Euler- 
like sequences are determined in this paper for the zero-pretwist case, as well as 
for the case of initially applied pretwist and the case of combined pretwist and 
deformation twist. In addition, it is shown that twist variables associated with the 
combined twist case can be related to those for which the pretwist is applied before 
deformation (examples are included). Using these relationships, the different forms 
of the transformation matrix can be made identical by performing a change of twist 
variable. 

The forms of the curvature vector associated with the different sequences and 
pretwist treatments can also be reduced to one form by performing a similar change of 
twist variable. Relations between the elastic twist and the twist variables associ- 
ated with the flap-lag-twist sequence, lag-flap-twist sequence, and the alternate 
approach are listed. 

By substituting expressions for one twist variable in terms of another, it is 
hoped that equations of motion developed using different twist treatments can be com- 
pared. Naturally, other differences in the derivation steps may also have to be 
dealt with in making any two sets of equations comparable. 
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TABLE 1.- SIX FORMS OF THE TRANSFORMATION MATRIX FOR THE NO-PRETWIST CASE 


Flap-lag-twist 


W x ) = 


-(ab cos^ + c simf^) 


-J 75 (1 - b 2 ) 1 / 2 cos* / 

) ' 1 (1 - b 2 ) 1 ' 2 


(a sin^j - be cos*^ 


(1 - b- 


(ab sin*, - c cos*,) . (a cos*. + be sin*,) 

; - -a - b 2 ) 1 ' 2 sin*. - 1 J -■ 

(l - b 2 ) 1/2 1 (1 - b 2 ) 1 ' 2 


Flap-twist-lag 


19A 


-[ab(cos 2 * - b 2 ) 1 / 2 + c sin* 2 ] [a sin*, - bc(cos 2 * - b 2 ) 1 / 2 , 

2 (eos 2 * 2 - b 2 ) 1 / 2 ^ 

2 1 - b ? 

[a(cos 2 * - b 2 ) 1 / 2 + be sin* ? \ 


1 - b z 

[ab sin* ? - c(cos 2 * - b 2 ) 1 / 2 ] 


1 - b* 


-sine}) 


1 - b* 


Twist-flap-lag 




-a(b cos* 3 + c sin* 3 ) (1 - b 2 )cos*- - be sin*. (1 - c 2 )sin* 3 - be cos* 


( 1 - (b cos* 3 + c sin* 3 ) 2 ] 1 2 [1 - (b cos* 3 + c sin*^) 2 ] 1 / 2 [1 - (b cos* 3 + c sin* 3 ) 2 ] 1 /' 

-a sin* 


(b sin* 3 - c cos* 3 ) 


3 -- 

3' ? 

a cos*. 


[1 - (b cos* 3 + c sin4> 3 ) 2 ] 1 ^ 2 (1 - (b cos* 3 + c sin*^ 2 ] 1 / 2 [1 - (b cos* 3 + c sin* 3 ) 2 ] */' 


La g-f lap-twist 


[&} 


-(b cos*. + ac sin* ) (a cos* 4 - be sin* 4 ) 


(1 - c 2 ) l/2 


(1 - c 2 ) 1 ' 2 


(b sin* 4 - ac cosf ) -(a sin* 4 + be cos* 4 ) 


(1 - c 2 ) 1 ' 2 


(1 - C z ) 


2xl/2 


(1 - c 2 ) 1 / 2 sin* 4 


(1 - C 2 ) 1 ^ 2 COS* 


Lag-twist-flap 



a 

b 

c 



-[b(cos 2 * 5 - c 2 ) 1 / 2 + ac sin* 5 

[a(cos 2 * 5 - c 2 ) 1 ' 2 - be sin* 5 ] 

sin* 5 


[$;] « 

1 - c 2 

1 - c 2 



[b sin* 5 - ac(cos 2 * 5 - c 2 ) 1 ^ 2 

-[a sin* s + bc(cos z * 5 - c 2 ) 1 ^ 2 ] 

(cos 2 * 3 - c 2 ) 1 /* 



1 " c 2 

1 - c 2 


Twist-lag-f lap 





a 

b 

c 

- 


-(b cos** + c sin**) 

a cos** 

a sin** 


[^»] = 

[1 - (b sin** - c cos**) 2 ] 1 / 2 

(1 - (b sin** - c cos**) 2 ] 1 ' 2 [1 

- (b sin** - e cos**) 2 ] 1 / 2 


a(b sin* 6 - c cos**) 

-(1 - b 2 )sin** - be cos** 

(1 - c 2 )cos* 6 + be sin** 


[1 - (b sin** - c cos**) 2 ] 1 / 2 

[1 - (b sin** - c cos**) 2 ] 1 / 2 [1 

- (b sin** - c cos**) 2 ] 1 / 2 
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TABLE 2.- RELATIONSHIPS BETWEEN TWIST VARIABLES FOR THE NO-PRETWIST CASE 
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TABLE 3.- FOUR FORMS OF THE TRANSFORMATION MATRIX FOR THE CASE OF INITIALLY APPLIED 

PRETWIST 


Flap-lag-twist 


[-a(b cos0 + c sin0)cos4 ()l {-a sinQ sin<lQ 1 {a cos8 sinct^ 

+ (b sinO - c cosO)sin<ftQ ] + ( (1 - b 2 )cos0 - be sin0 ]cos^ lJl } + [ (1 - c 2 )sin0 - be cosO ] cos<f>g x } 

[^ ] = | 1 - (b cosG + c sinG) 2 ] 1 ^ 2 [1 - (b cos0 + c sinO) 2 ] 1 ^ 2 [i - (b cosQ + c sinG) 2 ] 1 / 2 

{a cosO cos^Qj 


i (b cos0 + c sin0)sin^0 1 {-a sin0 


+* (b sin6 - c cos0)cos0q 1 J - [(1 - b 2 )cos0 - be sinO]sin^Q 1 } - [(1 - c 2 )sin0 - be cosO ]cos<J>q 1 } 

{1 - (b cos 9 + c sinO) 2 ] 1 / 2 [1 - (b cos 0 + c sinG) 2 ] 1 / 2 [1 - (b cos 0 + c sinG) 2 ] 1 / 2 


[1 - (b cos 0 + c sinG) 2 ] 1 / 2 


Twist-flap-lag 


-a(b cosa 3 + c sina 3 ) 


\@ ] = U - (b cosa 3 + c sina 3 ) 2 ] ] 


(1 - b 2 )cosot 3 - be sina 3 
[1 - (b cos*'t 3 + c sinu 3 ) 2 ] 1 ^ 2 


(1 - c z )sina 3 - be cosa 3 
[1 - (b cosa 3 + c sinou) 2 ] 1 / 2 


(b sina 3 - c cosa 3 ) 

[1 - (b cosa 3 + c sina 3 )‘ L ] 1 / 2 


[1 - (b cosa 3 + c sina^ 2 ] 1 / 2 


[1 - (b cosa 3 + c sinag) 2 ] 1 / 2 


Lag-f 1 ap-twist 



[-(b cosO + c sinO)cos<|) 0 4 
F a (b sin0 - c cosO) sin4 >(.) 4 ] 

{a cosO cos<J>q 4 
- [ (1 - b 2 )sinO 

+ be cosO]sin 4 ) 0 4 } 

{a sinG cos(j>Q 4 
+ [ (1 - c 2 )cos 0 

+ be sinO ] sin<J)g 4 } 

CD ' 

-r 

II 

[1 - (b sinG - c cosO ) 2 ] 1 / 2 

[1 - (b sinO 

- c cosO ) 2 ] 1 / 2 

[1 - (b sinG 

- c cosO ) 2 ] 1 / 2 


[ (b cos 0 + c sin 0 )sin 4 »Q 4 
+ a(b sin 0 - c cos 0 )cos<j> 0 4 ] 

{-a cosO sin<j> 0 4 
- ((1 - b 2 )sin 0 

+ be cosO ] cos<J>q 4 } 

{-a sinG sin<$> 0 4 
+ [ 0 - c 2 )cos 0 

+ be sinG jeos^Q^} 


[1 - (b sinG - c cosO ) 2 ] 1 / 2 

[1 - (b sinO 

- c cosO ) 2 ] 1 / 

[ 1 — (b sinG 

- C cose ) 2 ] 1 / 2 


! Twist-lag-flap 


-(b cosa 6 + c sina 6 ) 

(&.J » [1 - (b sina 6 - c coso^) 2 ] 1 / 2 


a(b sina g - c cosa 6 ) 

[1 - (b sina £ - c cosag) 2 ] 1 / 2 


[1 - (b sina 6 - c cosa^) 2 ] 1 / 2 
-(1 - b‘')sina g - be- cosa^ 

[1 - (b sina 6 - c cosa 6 ) z ] x / 2 


[1 - (b sina 6 - c cosa 6 


(1 - c z )cosa 6 + be sina 6 
[1 - (b sina 6 - c cosag) 2 ] 1 / 2 
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TABLE 4.- RELATIONSHIPS BETWEEN TWIST VARIABLES FOR THE CASE OF INITIALLY 

APPLIED PRETWIST 


r 


t.an<j> 


02 


tan<}> 


(1 - B 2 ) l/2 tan(j ) 01 
(1 + B 2 tan 2 <J > 01 ) l/2 

01 


To 

<P 


0(e 2 ) 


(1 - B 2 )tan<(> 


03 


A + BC taring 
A taring i - BC 


02 01 


*63 ’ *61 


| tan< ^04 A + BC tan(j) 


01 


tan<}> 


A tan<j > 01 - BC 


K 0 s 


tan<j> 


06 


r 

tan<|> g g 


[(1 - B 2 ) - B 2 C 2 + 2 ABC tan«fr 01 + C 2 tan 2 <f> 01 ] x / 2 

A tan$ 0 BC 
(1 - B 2 ) 

(1 - C 2 ) l/2 taring 


<h s <b - BC 
y 04 V 0l 


*05 S *01 " BC 


*06 “ *01 " BC 


(1 + C 2 tan 2 ())g 4 ) 1 


71 


(1 - C 2 )tan<f> 


tan<t> 


tan<p 


04 


06 


03 


A - BC tan<J > 04 

(1 - B 2 )tan<j> 


*05 * *04 


*06 = *04 


02 


BC tan <{> e2 + A[1 - B 2 (1 + tan 2 (j)g 2 ) ] 1 / 2 
(1 - C 2 )tan<t>, 


tan *e$ 

.where 


P 05 


-BC tan<f > 5 + A[1 - C 2 (1 + tan 2 $ 0 S ) ] x / 2 


A = a 

B = b cos0 + c sin0 
C = -b sin0 + c cos0 
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TABLE 5.- SIX FORMS OF THE TRANSFORMATION MATRIX TO 0(e) FOR THE COMBINED-TWISTS 

CASE 


Flap-lag- twist 

r 1 - b 2 /2 - c z /2 

-(b cosO + c sinO) 

+ 0 ai (b sinO - c cost)) 

(b sinO - c cos0) 

+ <J> (b cosO + c sinO) 


Flap-twist-lag 


[ff 1 

1 '■'oil* 


1 - b 2 /2 - c 2 /2 


\<y } 

1 a 3 J 




Lag-twist-flap 




1 - b-/2 - cr/2 


Twist-lag-flap 


f.sr 


1 - b 2 /2 - c 2 /2 


(4> - bc)cosO 


(1 - c" /2 - ♦* i /2)sinO 


-(b cosO + c sinO) 

+ 4> a2 (b sinO - c cosO) 

(b sinO - c cosO) 

L+ 4> u2 (b cosO + c sinO) 
Twist- flap- lag 

1 - b 2 /2 - c 2 /2 

-(b cosO + c sinP) 

■f 4>ot 3 ( b sinO - c cosO) 

(b sinO - c cosO) 

+ 4> a3 (b cosO + c sinO) 

Lag-flap-twist 

1 - b 2 12 - c 2 /2 


(1 - b 2 /2 - 4 > 2 1 /2)cos0 - 4> al sinO + 

-U\ tl - be) sinO 

-<1 - b 2 /2 - ♦ 2 1 /2)sin0 - « (U cosO + (l _ c 2 /2 _ ^ /2)co!lC 


(1 - b"j2 - <{) 2 2 /2)cosO ■ (4» w2 - be) cost* 

- (<J> a2 + b 2 /2 tanG)sinO + (1 - c 2 /2 + b"/2 - 4>^ 2 /2)sinO 


-O - f z /2)sin0 - 6 cos6 


-(4> ft , - be - b 2 /2 tanO)sinO 
+ (1 - c 2 /2 - * 2 /2)cosO 


-{1 - l/2(b sinO - c cosO) 2 - 4> z 3 /2]sinO 
- <J> a3 cosO 


+ [1 - c‘ + l/2(b cosO + c sinO)*' ~ <£ /21sinO 

(i 3 

+ (1. - 1/2 (b sinO - c cosO) 2 - 4^ /2]cosO 


-(b cosO + c sinU) 

+ * <b sinO - c cosO) 

(b sinO - c cosO) 


+ be) sinO 

+ (1 - b z /2 - #* % /2)cosU 

-(•!>,, n + be) cosO 
- (1 - b z /2 - .» z /2)«inO 


(1 - e 2 / 2 - (,V /2)siiv + o cosO 


(1 - e 7 /2 - ;• /2)cos - sine 


^(b cosO + c sinO) 

+ 4)^ 5 (b sinO - c cosO) 

(b sinO - c cosO) 

4>, t5 (b cosC + c sinO) 


-(4* l5 + be + e 2 /2 tan'OsintJ 
+ <1 - b z /2 - 4> 2 5 /2)cost) 

-(4* r - bc)cosO 


(1 - 40 /2) si n-- + »: oos 4 * 


(1 - c /2 - B /2)ros- 


-(b cosO + c sinO) 


(b sinO - c cosO) 

L + g (b cosU + c sinU) 


- (1 - b z /2 - c z /2 - 4« 2 5 /2)sin<> - (y l5 + c z /2 tan )sin- 


-* a6 sinO fl - 1/2 (b cos*; + c sin )’ - 4* 2 6 /2 ]sint: 

+ [l - 1/2 (b cost) + c sinO)" - 4 2 /2|cos0 + ■: , cosv 

(iS '’® 

~(0 ic - + bc)cosO [] - c 2 + 1/2 (b sin - - c cosO)*' - /21cos( 

•*® 1 6 

- [1 - b 2 + l/2(b sinO - c. cosO) 2 - 4 O ./2]sinO - (4- , - bc)sin* 
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TABLE 6. 


tana 2 


tana 3 


< tana 4 


tana g 


j^tanOg 

tana x 


tana 2 


tana. 

4 


tano 5 

where 


- RELATIONSHIPS BETWEEN COMBINED-TWIST VARIABLES 
(I - b 2 ) l/z tanc^ 

(1 + b 2 tan 2 a 1 ) 1 ^ 2 

(1 - b 2 )tana i 
a + be tana x 

a tana i - be 
a + be tana x 

a tano^ - be 

[(1 - b 2 ) - b 2 c 2 + 2abc tana x + c 2 tan 2 a 1 ] 1 ^ 2 
a tano^ - be 

a - b 2 > 

a tana g 

r (1 - b 2 ) - be tana 3 

a tana, 

_ 3 

[(1 - b 2 ) - 2bc tana 3 + b 2 tain 2 a 3 ] 1 ^ 2 

a tana 

6 

(1 - c 2 ) + be tana 6 

a tana, 

6 

[(1 - c 2 ) + 2bc tana 6 + c 2 tan 2 a g ] 1 ^ 2 

a i = 9 + *ie 

a 2 = 0 + *20 

a 3 = 6 + *30 

% = 6 + *40 
a 5 = 0 + <f> 50 

a e = 0 + *60 
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TABLE 8.- RELATING TWIST VARIABLES ASSOCIATED WITH DIFFERENT PRETWIST 
TREATMENTS FOR THE FLAP -LAG-TWIST AND LAG-FLAP-TWIST ROTATION 
SEQUENCES 


a tana 3 

tana = 

(1 - b 2 ) - be r.ana 3 

[a + be tan 0 + (1 - b 2 )tan 2 0 ]tan<b - [(1 - b 2 - a)tan 0 - be tan 2 0 ] 

:an<p - 

ai [(1 - b 2 ) - be tan 0 + a tan z 0 ] - [be + (1 - b 2 - a)tan 0 ]tan<j > a3 

[(1 - b 2 ) - be tan 0 + a tan 2 0 ]tan<f> - [(1 - b 2 - a)tan 0 - be tan 2 0 ] 

, • 01 
:an<}> = 

a ‘ [ (1 - b 2 ) - be tan 0 + a tan 2 0 ] + [(1 - b 2 - a)tan 0 - be tan 2 0 ]tan<}) 0 :i 

4>ai = 4 'q i + j (b 2 - c 2 )sin 0 cos 0 + be sin 2 0 to 0(e 2 ) 


a tana, 

6 

(1 - c 2 ) + be tana g 

[a - be tan0 + (1 - c 2 )tan 2 0]tan<j> ag - [(1 - c 2 - a)tan0 + be tan 2 0] 
[(1 - c 2 ) + be tan0 + a tan 2 0] - [-be + (1 - c 2 - a)tan0]tan<t> ag 

[ (1 - c 2 ) + be tan0 + a tan 2 0 ]tan<|> Q — [ (1 — c 2 — a)tan0 + be tan 2 0] 
04 

[ (1 - c 2 ) + be tan0 + a tan 2 0] + [(1 - c 2 - a)tan0 + be tan 2 0]tan4> e . 


ij> 0 4 - -j (b 2 - c 2 )sin 0 cos 0 - be sin 2 0 to 0 (e 2 ) 
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TABLE 9.- RELATIONSHIPS BETWEEN THE ELASTIC TWIST AND SEVERAL DEFORMATION-TWIST 

VARIABLES FOR THE NO-PRETWIST CASE 



28 



TABLE 10.- RELATIONSHIPS BETWEEN THE ELASTIC-TWIST AND SEVERAL DEFORMATION-TWIST 

VARIABLES FOR THE NONZERO-PRETWIST CASE 


(❖’ + e 


_ dx r i + a‘0* + (b cosO + c sinG)(b f sinO - c 1 cosO) + (b cosO + c sint)) x (c , b - b'c) 1 
^ X 3 [ 01 a [ 1 - (b cosO + c sinO) 2 ] J 

+ T X f , [ Vo' + (b cosO + c sin 1 *) (b ' sinO - c' cosO) + (b cosO + c sinU) 2 (c , b - b'cH l 
* Gl " e J 0 \ L a [ L - (b cosO + c sinO) 2 ] Jj 


? <+> e + J f(b cosO + c sinO)(-b sinO + c cosO)' + — (b“ + c ? )0‘Jdx 


to 0(e) 




a (=' *r?V b ')'‘" 

= 4> e + J be' dx to 0(e 2 ) 


— (4>I + G') = + 0’ + - 

x dx 04 a 


^4 = *e 


!x /lf , ntx dx T . a 2 0' + (-b sinO + c cos0)(b* cosO + c* sinO) + (-b sinO + c cosO) ? (c'b - b'c)l 

= ^7 L°" 7 a[l - (-b sine + c cosO) 2 ] J 

n , I" a 2 0 * + (-b sinO + c cosO)(b' cosO 4- c' sinO) + (-b sinO + c cosu) 2 (c'b - b'c) ~| 1 

L a [ 1 — (-b sinO + c cosO) 2 ] J / 

f f(-b sinO + c cos0)(b cosO + c sinO)' - y (b 2 + c 2 )0'ldx 

<o *- 

[*« * •' *•(“■' *rr?'')] 

j b’c dx to 0(c 2 ) 


t ( *e + 0 ') = ^ L^ + 6 


+ e- + 


■t> T e _ 


♦f 


f X / be ' - b 'c \ 


dx to 0 (r 2 ) 







ai 


; i = 1,2, 3. 4, 5, 6 


SUBSTITUTE FOR sin a,, cos a-, IN TERMS OF sinofj, cosaj WHERE 

tan aj = f j (a, b, c, tan a-.) ; i = 1 , 2, 3, 4, 5, 6 

j - 3, 6 


r‘J 

♦ 

SUBSTITUTE FOR sinaj, cosaj WHERE 
a j = 6 + 0<xj = e + ^0j - j = 3 ' 6 

T e\ 

♦ 

APPLY TRIGONOMETRIC IDENTITIES 

cosctj = cos 0 cos j - sin0 sin 00 j ; j = 3, 6 
sinaj = cos 0 sin + sin 0 cos 00 j 

4 

SUBSTITUTE FOR sin 4> e - cos0 0 j IN TERMS OF sin <t> ek , cos0 0k WHERE 

tan 00: = f (a, b, c, 0, tan 0 0k ) ; j = 3, 6 

k = 1,2, 3, 4, 5,6 


r 0k 


Figure 5.- A series of variable changes between transformation matrices 
associated with the combined twists and initial pretwist treatments. 
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Figure 7.- Bend-twist and twist-bend rotation sequences for the case of no pretwist. 
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